Geometry of random interactions 
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It is argued that spectral features of quantal systems with random interactions can be given a 
geometric interpretation. This conjecture is investigated in the context of two simple models: a 
system of randomly interacting d bosons and one of randomly interacting fermions in a j = | shell. 
In both examples the probability for a given state to become the ground state is shown to be related 
to a geometric property of a polygon or polyhedron which is entirely determined by particle number, 
shell size, and symmetry character of the states. Extensions to more general situations are discussed. 



Recent studies in the nuclear shell model 0, 0, 0> 3 an d 
the interacting boson model 0, ||J with random interac- 
tions have unveiled a high degree of order. In particu- 
lar, a marked statistical preference was found for ground 
states with J = + . In this Letter it is argued that 
spectral features of quantal systems with random inter- 
actions can be given a geometric interpretation which 
allows the computation of the probability for the quan- 
tum mechanical ground state to have a specific angular 
momentum, based on purely geometric considerations. 
Although these results are obtained in the context of a 
variety of simple models which do not cover the full com- 
plexity of random interactions, we believe them to be 
sufficiently general to conjecture the possibility of an en- 
tirely geometric analysis of the problem. 

Consider a system consisting of n interacting particles 
(bosons or fermions) carrying angular momentum j, inte- 
ger or half-integer. Eigenstates of a rotationally invariant 
Hamiltonian are characterized as \j n aJM) where J and 
M are the total angular momentum and its projection, 
and a is any other index needed for a complete labeling 
of the state. Although spectral properties of a Hamil- 
tonian with both one- and two-body interactions can be 
analyzed in the way explained below, we assume for sim- 
plicity that the one-body contribution is constant for all 
eigenenergies and that the energy spectrum is generated 
by two-body interactions only. Under this assumption its 
matrix elements can be written as 



(j n aJ\H\j n a'J) 



n(n — 1) 



(1) 



where M is omitted since energies do not depend on 
it. The quantities Gl are two-particle matrix elements, 
Gl = (j 2 L\H 2 \j 2 L), and completely specify the two- 
body interaction while c^ aJ are interaction-independent 
coefficients. They can be expressed in terms of coeffi- 
cients of fractional parentage (CFP) Q and, as such, are 
entirely determined by the symmetry character of the n- 



particle states. 

We begin by considering the special case when a ba- 
sis \j n aJM) can be found in which the expansion Q is 
diagonal. In this case the Hamiltonian matrix elements 
reduce to the energy eigenvalues 



n(n - 



L 
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with6^ Qj = (c^ qJ ) 2 . This is obviously a simplification of 
the more general problem Q but nevertheless a wide va- 
riety of simple model situations can be accommodated by 
it. For example, this property is valid for any interaction 
between identical fermions if j < | and remains so ap- 
proximately for larger j values; it is also exactly valid for 
p, d, or / bosons. We shall refer to this class of problems 
as diagonal. For a constant interaction, Gl = 1, all n- 
particle eigenstates are degenerate with energy |n(n — 1) 
and consequently the coefficients b^ aJ satisfy the prop- 
erties J2l bnaj = 1 an d < b^ a j < 1. Equation J3J) can 
thus be rewritten in terms of scaled energies as 
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for arbitrary L' . This shows that, in the case of N in- 
teraction matrix elements Gl, the energy of an arbitrary 
eigenstate can, up to a constant scale and shift, be repre- 
sented as a point in a vector space spanned by m = N — 1 
differences of matrix elements. Note that the position 
of these states is fixed by b^ aJ and hence interaction- 
independent. Furthermore, since < b^ aJ < 1, all states 
are confined to a compact region of this space with the 
size of one unit in each direction. For independent vari- 
ables Gl with covariance matrix (GlGl 1 ) — Sll'i states 
are represented in an orthogonal basis. The differences 
in (PI are not independent but have a covariance matrix 
of the form 1 + Sll 1 ', this leads to a representation in a 
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m-simplex basis (i.e., an equilateral triangle in m = 2 
dimensions, a regular tetrahedron for m = 3,. . . ). 

The following procedure can now be proposed to de- 
termine the probability P na j for a specific state naJ to 
become the ground state. First, construct all points cor- 
responding to the energies e na j. Next, build from them 
the largest possible convex polytope (i.e., convex poly- 
hedron in m dimensions). All points (i.e. states) inside 
this polytope can never be the ground state for whatever 
choice of matrix elements Gl and thus have P na j = 0. 
Finally, the probability P na j of any other state at a ver- 
tex of the polytope is a function of some geometric prop- 
erty at that vertex. 

This general procedure can be illustrated with some 
examples. Consider first a system of j = 2 d bosons. 
In this case there are three interaction matrix elements, 
G , G 2 , and G 4 with G L = (d 2 L\H 2 \d 2 L); thus, N = 3 
and the problem can be represented in a plane. Because 
of the U(5) D SO(5) D SO(3) algebraic structure, an 
analytic solution is known for n interacting d bosons with 
eigenenergies H,E1,E3 



2n(n - 2) + 2i>(v + 3) - J(J + 1) 



n(n + 3) 



7n(n — 1) 



5n(n — 1) 



(G - G 4 ) 



(G2 — G4) 
G 4 , (4) 



where v is the seniority quantum number which counts 
the number of d bosons not in pairs coupled to angular 
momentum zero. Any energy e nu j can be represented 
as a point inside an equilateral triangle with vertices Go, 
G2, and G4 of which the position is determined by the 
appropriate values on the edges Go— G4 and G2—G4. Ex- 
amples for several boson numbers n are shown in Fig. ^ 
For n — 2 there are three states with L = 0, 2, 4 and ener- 
gies e = Go, G2, G4; clearly, they have equal probability 
of being the ground state. As n increases, more states ap- 
pear in the triangle. The majority of states, shown as the 
smaller dots, are inside the convex polygon indicated in 
grey and can never be the ground state for whatever the 
choice of Gl. If we translate or rotate the convex poly- 
gon inside the triangle, its properties should not change 
since the points Go, G2, and G4 are equivalent and since 
the distribution depends only on G 2 . Thus the prob- 
ability for a point to be the ground state can only be 
related to the angle subtended at that vertex. The rela- 
tion can be formally derived but also inferred from a few 
simple examples. If the polygon is an equilateral triangle, 
square, regular pentagon,. . . each vertex is equally proba- 
ble with probability ~, |, . . One deduces the relation 
(see the m = 2 polygons in Table |lj 



p(2) = i _ fl 
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between the angle 9 V at the vertex v of the polygon and 

(2) 

the probability P„ ; for the state associated with that 




FIG. 1: Polygons corresponding to several systems of interact- 
ing d bosons. All states are represented by a dot. The smaller 
dots are inside or on the edge of a convex polygon (with the 
larger dots as vertices and indicated in grey) and are never 
the ground state. The probability of the vertex states to be 
the ground state is related to the angle at the vertex. 



TABLE I: The angle or the angle sum X^/gi, an d the ver- 
tex probability P v of regular polygons (m = 2) and polyhe- 
drons (m = 3). 
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vertex to be the ground state. Table HT| compares prob- 
abilities calculated with the analytic relation (|SJ with 
those obtained from numerical tests for several systems 
of randomly interacting d bosons with boson numbers 
n = 5, 6, 10, 18. The numerical probabilities are obtained 
from 20000 runs with random interaction parameters. 
They agree with the analytic result JSJ. 

A second example concerns a system of four j = 
I fermions which was discussed recently by Zhao and 
Arima Q. In this case there are four interaction ma- 
trix elements, G , G 2 , G4, G 6 with Gl = {j 2 L\H 2 \j 2 L}, 
and this leads to a problem that can be represented in 
three-dimensional space. Any state in the j = | shell 
can be labeled with particle number n, seniority v, and 
total angular momentum J with analytically known ex- 
pansion coefficients b^ vJ 0. For n = 4 there are eight 

nuj 



different states and the corresponding coefficients b. 
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TABLE II: Probabilities P nv j (in %) of some d-boson states 
obtained analytically with formula and numerically with 
Gaussian parameters. 
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10 


0(0) 


20.87 
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0.50 


0.41 




2(5) 


36.19 


36.59 




0(10) 


37.56 


37.96 




10(5) 


39.63 
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38.37 




12(6) 


39.63 


39.49 
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41.75 



are given in Ref. |3j. These eight states can be repre- 
sented in three-dimensional space spanned by the three 
differences Go — Gq, Gi — Gq, and G4 — Gq, and six of 
them define a convex polyhedron (see Fig. 0). The two 
remaining points [corresponding to the state J{v) = 4(2) 
and 5(4)] are inside this polyhedron and are never the 

ground state. The relation between the geometry of the 

(3) 

polyhedron and the probability Py ' of each vertex state 
to be the ground state can again be inferred from a few 
simple examples. The relevant quantity in this case is 
Ylf3v wnere the sum is over all faces that contain 
the vertex v and 9 v f is the angle at vertex v in face /. 
A few examples with regular polyhedrons (see the m = 3 
polyhedrons in Table [IJ demonstrate that the relation is 

f3v 

Table IIIII compares the probabilities for the different 
states to become the ground state as calculated in var- 
ious approaches. The second column gives the analytic 
results obtained from © while the third column lists nu- 
merical results obtained from 20000 runs with random 
matrix elements with a Gaussian distribution. As our 
code does not distinguish between states with the same 
angular momentum J but different seniority v, only the 
summed probability for each J is given. The last column 
shows the results of Zhao and Arima 3] who calculate 
the probability as a multiple integral. 

These notions can be generalized in several ways. The 
first is towards energies that depend on a set of con- 
tinuous variables {ti,...,t q } as follows [compare with 
Eq. ©]: 

E(t u ...,t q )=J2 b L(tu---,t g )G L . (7) 

L 

The analogous problem now consists in the determina- 
tion of the probability density dP{t\, . . . , t q ) for obtain- 
ing the lowest energy at {t\, . . . ,t q } with random inter- 
actions Gl- We assume by way of example that the 




FIG. 2: The polyhedron for a system of n = 4 j = i fermions. 
The upper part indicates the position of the polyhedron in the 
tetrahedral coordinate system. The left part shows the same 
(enlarged) polyhedron with vertices that correspond to the 
states J{v) that can become the ground state with a proba- 
bility given by the exterior angle. Two additional states with 
J{u) — 4(2) and 5(4) lie inside or on the face of the polyhedron 
and are never the ground state. For representation purposes 
the polyhedral face 0(0)-2(4)-8(2) has been removed in the 
left part. 



TABLE III: Probabilities P mJ (in %) of n = 4 j = | states 
obtained with the analytic formula JSJ , from a numerical cal- 
culation, and from the integral representation of Zhao and 
Arima Q. 
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number of variables q is one less than the number N 
of interactions Gl, q = N — 1. In that case Eq. (JJJ rep- 
resents a g-dimensional hypersurface Tfl embedded in a 
(q + l)-dimensional Euclidean space E 9+1 (the metric fol- 
lows from the covariance matrix (GlGl 1 ) = 6ll')- Let 
us suppose that £ 9 is the closed orientable manifold. It 
can then be shown that the probability density is given 
by Gauss' spherical map [Jjj £ 9 — > S q where S q is a q- 
dimensional hypersphere of a unit radius. If the degree 
of the spherical map is one (as it is for closed convex 
surfaces), the probability density reads 

dP{h,...,t q ) = —K q (t 1 ,...,t q )dv, (8) 

where K q is the Gaussian curvature of £ 9 , S q is the vol- 
ume of the unit hypersphere S q and dv denotes an in- 
finitesimal element of £ 9 . In the simplest example of 
one parameter t\ = t and two interactions G\ and G2, 
the energy is parametrized as a curve in a plane. It can 
then be shown that the probability density is given by 
dP(t) = (27r) _1 Ki(t)ds, where ds is the infinitesimal arc 
length. In fact, this formula is also valid for piecewise 
curves and precisely leads to the result (J5J for a polygon. 
The validity of the result can be checked by comparing 
the probability obtained by integration of {8J over a part 
of Yj q to the numerically calculated one. As an exam- 
ple we discuss a two-dimensional ellipsoid in E 3 with one 
semi-axis c different from the others a. The probability 
associated with the part of the surface with spherical co- 
ordinates (9, 4>) satisfying < 6 < Qq and < <fr < <fio is 
given by 

p(e ,<p ) = ^(i- , aWo V (9) 

4tt y y/c 2 + (a 2 - c 2 )(cosc/>o) 2 / 

This expression has been compared with the numerically 
calculated probability; the difference is close to zero. We 
have analyzed likewise the case of a three-dimensional 
hyperellipsoid in E 4 , showing that our approach can be 
generalized to higher dimensions. These results also open 
up the possibility for an extension to higher-dimensional 
polytopes, by replacing the right-hand side of (Eq. |8J) 
with an appropriate characteristic for each vertex of the 
polytope. Indeed, it can be shown that the probabilities 
Py 771 ^ in © and © are related to the exterior angle at 
vertex v [l2.ll3| of either a convex polygon (m = 2) or a 
convex polyhedron (m = 3). 

We believe that, although derived for a restricted form 
of interaction Hamiltonians, these results suggest that 
generic rt-body quantum systems, interacting through 
two-body forces, can be associated with a geometrical 
shape defined in terms of CFPs or generalized coupling 
coefficients. Geometry thus arises as a consequence of 
strong correlations implicit in such systems and is in- 
dependent of particular two-body interactions. Random 



tests can be understood in this context as sampling ex- 
periments on this geometry. In this Letter we have shown 
that geometric aspects of n-body systems determine some 
of their essential characteristics. In particular, for diago- 
nal Hamiltonians surface curvature defines probability to 
be the ground state. Other correlations could also be re- 
lated to geometric features. These results generalize and 
put onto a firm basis the previous work which hinted at a 
purely geometric interpretation of randomly interacting 
boson systems 0, , as well as provide an explanation 
for the method of Zhao and collaborators |3|. In fact, 
in the latter reference, the authors have advanced some 
qualitative reasons for certain states to dominate and 
later provided an approximate procedure (which is not 
always accurate to estimate the ground state prob- 
abilities, although no reason was offered for its success. 
Our study, at least for the case where the Hamiltonian 
is diagonal, clearly shows that the procedure of Zhao et 
al is equivalent to a projection of the considered above 
polyhedra on the axes defined by the two-body matrix 
elements, which tend to correlate well with the angles we 
introduce. This connection will be elaborated in detail 
elsewhere |l5j |. Further work is required to fully explore 
the geometry and its consequences for our understanding 
of n-body dynamics. 
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